Abstract. The narrow- line Seyfert 1 galaxy Mrk 766 was 
observed for 60 ks with the XMM-Newton observatory. 
The source shows a complex X-ray spectrum. The 2-10 
keV spectrum can be adequately represented by a power 
law and broad Fe Ka emission. Between 0.7 and 2 keV the 
spectrum is harder and exhibits a flux deficit with respect 
to the extrapolated medium energy slope. Below 0.7 keV, 
however, there is a strong excess of emission. The RGS 
spectrum shows an edge-like feature at 0.7 keV; the en- 
' ergy of this feature is inconsistent with that expected for 
. an OVII edge from a warm absorber. Mrk 766 varies by a 
' factor of 2 in overall count rate in the EPIC and RGS 
instruments on a timescale of a few thousand seconds, 
^ ■ while no significant flux changes are observed in the ul- 
O , traviolet with the OM. The X-ray variability is spectrally 
' dependent with the largest amplitude variability occur- 
, ring in the 0.4-2 keV band. The spectral variability can 
T— I ' be explained by a change in flux and slope of the medium 
■ energy continuum emission, superimposed on a less vari- 
able (or constant) low energy emission component. 
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1. Introduction 

X-ray emission currently provides the best tracer of condi- 
tions close to the supermassive black hole central engines 
of AGN. However the origin of AGN X-ray emission re- 
mains poorly understood. Between 1 and 20 keV AGN 
generally have a power law spectrum modified by the ef- 
fects of reflection and absorption (Pounds et al. 1990 ); be- 
low 1 keV many AGN also show some 'soft excess' emission 
(Walter & Fink 1993). The majority of recent models for 
AGN X~ray emission explain the soft excess as emission 
from an optically thick accretion disk, and the power law 
component as Compton upscattering of soft photons by a 
hot corona. 



Mrk 766 is a nearby {z ^ 0.0129), soft X-ray bright 
(0.5 - 2 keV flux of ^ 10~^^ erg s~^ cm ~^) a nd vari- 
able (factor 2 in ~ 1000s, Leighly et al. 1996) narrow 
line Seyfert 1 galaxy observed through a relatively small 
Galactic column of ~ 1.8 x 10^*^ cm^^. It is therefore an 
extremely good candidate for studying the dynamical and 
radiative processes which give rise to the X-ray emission 
in AGN. In this letter we present the spectrum and spec- 
tral variability of Mrk 766 as observed by the combination 
of instruments on board XMM-Newton during the perfor- 
mance veriflcation (PV) observation. 



Send offprint requests to: M.J. Page (mjp@mssl.ucl.ac.uk) 

* Based on observations obtained with XMM-Newton, an 
ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and the USA (NASA) 



2. The XMM-Newton observation and data 
reduction 

The observation was performed on 20th May 2000 during 
XMM-Newton revolution 0082 and lasted ^60 ks. All in- 
struments sufl[ered some data gaps close to the end of the 
observation. The data were processed using the standard 
analysis software (SAS) pipeline. For EPIC MOSl and 
PN analysis we use the most up to date (October 2000) 
on-axis response matrices available. The RGS data were 
analysed using response matrices generated by the SAS. 

Due to high particle background near the beginning 
of the observation, the EPIC PN and MOS cameras were 
closed for the first 20 ks and 30 ks respectively. The EPIC 
MOS cameras were operated in full frame imaging mode. 
The observed countrate of Mrk 766 is about five times 
larger than the recommended limit for this imaging mode, 
and therefore the MOS data are piled up. This is expected 
to result in a hardening of the MOS spectrum equivalent 
to a change in power law energy index of Aa ~ 0.05 (fig- 
ure 28 in the XMM-Newton User Handbook). Apart from 
this spectral hardening, we do not expect photon pile- 
up to seriously degrade the spectra. MOS spectra used 
in this paper were constructed using pattern 0-12 events. 
Source data were taken from a forty arcsecond radius re- 
gion around the source, and background data were taken 
from four identical circular regions close to the source po- 
sition. MOS Spectra were grouped to a minimum of 150 
counts per group. 

The EPIC PN camera was used in small window mode 
and is not significantly piled-up. Source data were ex- 
tracted from a 35 arcsecond radius region and background 
data were taken from three identical circular regions. Only 
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Table 1. Power law {Fe = K E'") fits to the EPIC 
spectra in the 2-5 and 7-10 keV energy ranges. 



Instrument 


a 


K* 




Ml 


0.98 ± 0.03 


4.66 ±0.16 


66.5/76 


M2 


0.96 ± 0.03 


4.30 ±0.16 


87.4/72 


FN 


1.09 ±0.02 


6.11 ±0.15 


71.8/71 


FN high 


1.14 ±0.03 


7.36 ±0.23 


12.7/18 


FN low 


1.05 ±0.03 


5.21 ± 0.18 


19.0/18 


* 10"=* keV 


cm~^ s~^ keV~ 


1 





single pixel (pattern 0) events were used to construct PN 
spectra, to match the available response matrix. In or- 
der to overcome a problem in current versions of the SAS, 
which produces periodic dips every '^25 eV in PN spectra, 
we have manually grouped PN spectra to groups which are 
> 25 eV. Less than 1% of the counts in the source regions 
in the MOS and PN cameras are due to the background. 

Both RGS cameras were operated in the standard spec- 
troscopy±Q mode for the duration of the observation. 
Source data were extracted using the standard spatial and 
order selections, and background data were taken from two 
identical sized extraction regions offset from the source in 
the cross dispersion direction. 87% of the counts in the 
source regions come from Mrk 766. 

The OM took a sequence of 5 exposures through the 
UVWl filter, 10 exposures through the UVM2 filter, and 
3 exposures through the UVW2 filter. 
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3. The spectra 

The EPIC PN, and MOS spectra of Mrk 766 are shown 
in Fig. and were analysed using the xspec software 
package. 



Fig. 1. a: spectra from MOSl (open squares), M0S2 
(open triangles), PN (filled dots) and power laws con- 
volved with the instrument response functions (stepped 
lines). b,c,d: ratios of MOSl, M0S2, and PN spectra to 
the power law model. 



3.1. 2-10 keV spectrum 

We began studying the spectra by fitting a power law 
model {Fe oc S"") with a Galactic Nh of 1.8 xlO^o cm'^ 
to the EPIC spectra in the energy ranges 2-5 and 7-10 
keV where the spectrum should be relatively unaffected 
by warm absorbers, Fe line emission or a soft excess. Table 
|l| shows the results of these fits. The data in these energy 
ranges are indeed fit well by a power law model. The best 
fit slopes for the MOS data are harder than those for the 
PN data by Aa '^0.1; this is probably due to pile-up in 
the MOS cameras (see Section |2|). The data and model 
(including the 0.2-2 and 5-7 keV regions not used in the 
fit) are shown in Fig. |^a. The ratio of the data to the 
model is shown for the three EPIC cameras in Fig. |l|b,c, 
and d. 

A prominent feature, consistent with a broad Fe Ka 
emission line centred at 6.4 keV, is present in the spectra 
from all 3 EPIC cameras. An emission line component is 
strongly required by the data: fitting the whole 2-10 keV 



PN spectrum with only a power law and fixed Galactic Nh 
results in a x^/^^ of 127/92, unacceptable at > 99% confi- 
dence. The presence of Fe Ka emission in Mrk 766 has so 



far been controversial: Leighy et al. (1996) found evidence 
for a narrow Fe Ka line, Nandra 1997 found evidence for 



a broad Fe Ka line (both using ASCA data), while Matt 
et al. 2000 found no strong evidence for an Fe Ka line at 
all in their BeppoSAX data. The results of fits to the 2-10 
keV spectrum with a power law continuum and Gaussian 
emission line model are given in Table ^, and residuals to 
these fits are shown in Fig. |2[ The Fe Ka line in Seyfert 
1 galaxies is widely thought to arise from fluorescence in 
near neutral material in an accretion disk around a black 
hole. An emission line produced by this process typically 
has an asymmetric profile with significant broadening to 
the red. The residuals to the gaussian fine fits (Fig. |2|) 
show an excess of emission between 5 and 5.5 keV, sug- 
gesting that the emission line in Mrk 766 has such a red 
wing. We have therefore investigated accretion disk emis- 
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Table 2. Fe Ka line parameters when the 2-10 keV spec- 
trum is fit with a power law -I- Gaussian emission line. 



Instrument 


Energy* 


a* 


EW* 




Model 


Ml 


6.5 ±0.2 


44+0-20 


r, oe+0.20 


78.2/88 


Energy" 


M2 


R 8+0-2 


u.oy_Q 21 




110/83 


EW 


PN 


6.5 ±0.1 




iQ+o-o-i 


90.4/89 


Inclination 

b 



Table 3. Relativistic accretion disk fits to the Fe Ka line 
in the PN data 

Schwarzschild Kerr 



* Rest frame in units of keV. EW is equivalent width. 




567567567 
Energy (keV) Energy (keV) Energy (keV) 

Fig. 2. Residuals of the gaussian fits to the Fe Ka line 
(Table |). Note the excess emission at 5-5.5 keV indicating 
that the emission line has a red wing. 



sion line models. Table gives the results of model fits 
to the Fe Ka line in the 2-10 keV PN data, which have 
the highest signal to noise. This time we have used emis- 
sion line profiles appropriate for accretion disks around a 
maximally rotating Kerr black hole (Laor 1991 , the laor 
mode l in x spec) and a Schwarzschild black hole (Fabian 
et al. 198£ , the diskline model in xspec). Both models fit 
the data well, and therefore the data are consistent with 
the Fe Ka line originating in an accretion disk. There is 
little to distinguish between the Kerr and Schwarzschild 
models in terms of ^ jv. 

Matt et al. (2000) found that fits to their BeppoSAX 
spectra of Mrk 766 were significantly improved with the 
addition of an edge at 7.55 keV, or an ionized refiection 
component (which produces an edge at a similar energy) . 
We have repeated fits to our 2-10 keV spectra with power 
law and emission line models, adding either an edge at 
~ 7.5 keV, or an ionized reflector (the pexriv model in 
xspec). For all three EPIC spectra, and for both addi- 
tional components, improvements to the are less than 
2(7 significant according to the F-test, and we therefore 
conclude that the spectra do not require an ionized reflec- 
tion component. 



3.2. 0.2-2 keV spectrum 

The EPIC data shown in Fig. ^ show that between 0.7 
and 2 keV the spectrum is harder than the power law 
extrapolation. Below 0.7 keV a strong excess is seen. Fig. 



6.5 ±0.1 
0.29 ±0.07 

34° ± 7° 

^■0-0.5 

86.3/88 



6.6 ±0.1 
0.37 ±0.11 

30° ± 6° 

Q 0+3.5 
'-'■^-0.4 

85.1/88 



" Rest frame in units of keV. EW is equivalent width, 
power law dependence of emissivity on radius: 
emissivity oc i?"**. Inner disk radius fixed at 6 Rg 



^ shows the spectrum of Mrk 766 recorded by the two 
RGS, in bins of 4 channels to improve signal to noise. 
The spectrum shows no strong narrow emission lines. This 
rules out a significant nuclear starburst contribution to 
the soft X-ray emission, which was suggested by Mathur 
( 2000| ). The most notable feature in the spectrum is a 
discontinuity at 0.7 keV. This feature corresponds to the 
highest energy of the soft excess emission seen in the EPIC 
spectra. A similar feature was previously identified with 
the K-edge of OVII from a warm absorber in ASCA and 
Beppo SAX data by Leighly et al. ( |1996D and Matt et al. 
( 2000D respectively; the presence of an OVIII K-edge was 
also inferred from both these datasets. At the redshift of 
Mrk 766 (0.0129) the OVII edge would be found at 0.730 
keV and the OVIII edge at 0.860 keV; these energies are 
labelled in Fig. ||. Edges are not seen in the RGS data at 
either of these energies. If the feature at 0.7 keV is an OVII 
K-edge, then the absorbing material must be redshifted by 
more than 10000 km s~^ which is inconsistent with current 
physical scenarios for the warm absorbe r gas , eg that it 
is in an outflowing wind (Kaastra et al. 2000 ) or that it 
is associated with the broad and/or narrow line regions 
(Otani et al. 1996 ). If the deficit of photons between 0.7 
and 2.0 keV is not produced by an ionized absorber, the 
continuum itself must deviate from a power law shape, 
and have a harder spectrum in this energy range than at 
higher energies. 

Detailed modelling of the soft X-ray spectrum of Mrk 
766 in this energy range is outside the scope of this let- 
ter. For a more detailed examination of the discrepancy 
between the Mrk 766 RGS data and warm absorber mod- 
els, see Branduardi-Raymont et al. ( |200l| ) who model the 
features in the RGS spectrum as emission lines from a 
relativistic accretion disk. 



4. Lightcurves 

RGS, EPIC PN and OM hghtcurves of Mrk 766 are shown 
in Fig. ^. The EPIC MOS lightcurves (not shown) are sim- 
ilar to the PN lightcurves but cover a shorter timespan and 
have lower signal to noise. The RGS and PN lightcurves 
are constructed with 300s bins. The PN lightcurves have 
been split into four different energy bands (0.1-0.4 keV, 
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Fig. 3. Spectra of Mrk 766 recorded by RGSl (top panel) 
and RGS2 (bottom panel). The lower line in each panel 
shows the size of the 1 a statistical uncertainty. The po- 
sitions of the expected warm absorber OVII and OVIII 
K-edges are labelled. 



0.4-2.0 keV, 2.0-6.0 keV, and > 6.0 kcV) and are shown 
in panels b-e of Fig. ^. Hardness ratios between the 0.1- 
0.4 keV, 0.4-2.0 keV and 2.0-6.0 keV energy bands are 
shown in panels g and h of Fig. ^. The low energy (0.4-2.0 
keV)/(0.1-0.4 keV) hardness ratio varies such that Mrk 
766 is harder when it is brighter. This is suggestive of a 
scenario in which the variability of a broad band emission 
component is diluted at the softest energies because of the 
superposition of a relatively constant soft excess compo- 
nent. Although the correspondence is not exact, for most 
of the observation the (2.0-6.0 keV)/(0.4-2.0 keV) hard- 
ness ratio appears to show the opposite behaviour to the 
lower energy hardness ratio: the 0.4-6.0 keV spectrum is 
harder when the 0.1-2.0 keV spectrum is softer and vice 
versa. 

We have quantified the variability in each of the FN 
lightcurves using the normalised variability amplitude 
(hereafter NVA) which we define as: 
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where a is the statistical uncertainty on the data points Xi . 
1 a statistical uncertainties on the NVAs were estimated 
by bootstrap resampling of the lightcurves. The NVAs of 
the 0.1-0.4 keV, 0.4-2.0 keV, 2.0-6.0 keV and >6.0 keV 
lightcurves are 0.15 ± 0.01, 0.19 ± 0.01, 0.17 ± 0.01 and 
0.12 ±0.02 respectively. The 0.4-2.0 keV lightcurve shows 



Fig. 4. Timeseries of Mrk 766 in the RGS (a), EPIC PN 
(b,c,d,e) and OM through the UVM2 filter (f). The y 
axis units for panels a-f are count s~^. Panels g and h 
show hardness ratios constructed from the PN lightcurves, 
demonstrating that the source undergoes significant spec- 
tral variability during the observation. 



more variability than the low energy (0.4-2.0) and high 
energy (>6.0 keV) fightcurves. 

The longest OM sequence of exposures (UVM2) is 
shown in Fig. In the ultraviolet there is little evidence 
for variability: a constant is an acceptable fit to the UVM2 
lightcurve [x^ /v of 11.1/9) and the data are all within 7% 
of the mean countrate. Similarly, no variability is found in 
the shorter UVWl and UVW2 sequences (not shown). 
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5. Spectral variability 

To examine how the spectrum of Mrk 766 changes with 
flux, we have divided the observation into two parts corre- 
sponding to "high flux" and "low flux" . These time periods 
are indicated in Fig. ^a. The EPIC PN spectra from the 
high and low periods are shown in Fig. after division 
by a model power law with energy index a = 1.09 (the 
best fit slope for the 2-10 keV PN spectrum of the whole 
observation) and fixed Galactic Nh- In order to improve 
signal to noise, the PN high and low flux spectra have been 
grouped more heavily than the PN spectrum for the whole 
observation. The main difference between the two spectra 
is that the overall continuum is higher in the higher flux 
spectrum. 

The ratio of the high to low flux spectra is shown in 
Fig. ||b. This shows that between 1 and 10 keV the high 
flux spectrum is softer than the low flux spectrum, while 
below 1 keV the high flux spectrum is harder than the 
low flux spectrum. We have performed power law fits to 
the PN data in the 2-5 and 7-10 keV energy range for the 
high and low flux spectra separately, and the results of 
these fits are given in Table |l|. The best-fit power law to 
the high flux spectrum is softer than the best fit to the 
low flux spectrum by Aa = 0.1, significant at 2a. Fig. ^ 
shows the ratio of the data to these best fit power law 
models. The two ratios are indistinguishable except be- 
low 0.7 keV, where the higher flux spectrum has a smaller 
fractional contribution from the soft excess emission. The 
overall spectral variability is therefore well characterised 
by an X-ray continuum which softens as it increases in 
flux combined with a less variable or constant compo- 
nent at low energy. This scenario would naturally explain 
the behaviour of the hardness ratios and why the 0.4-2.0 
keV lightcurve is the most variable of the PN lightcurves 
(Section ^). A variable warm absorber is not required to 
account for the spectral variability (as was suggested by 



Matt et al. 2000), because the deficit of photons between 
0.7 and 2 keV relative to the extrapolated higher energy 
power law spectrum is the same in the high and low flux 
spectra. 

The spectral variability between 0.3 and 2 keV can be 
examined at higher resolution using the RGS data. The 
high and low flux RGSl spectra of Mrk 766, in 16 channel 
bins, are shown in Fig. ^,b. Fig. ^ shows the difference 
between these two spectra, and Fig. ^ shows the ratio of 
the high flux to low flux spectra. As is found for the PN 
spectra, the ratio of the high flux to low flux RGS spectra 
is smaller close to 0.3 keV than above 1 keV. The 0.7 
keV edge-like feature is present in the low flux spectrum 
and the high flux spectrum. It also appears to be in the 
difference spectrum (though not the ratio spectrum) and 
hence if this feature is produced by a relativistic e missio n 
line, as proposed by Branduardi-Raymont et al. (2001), 
the line is varying with the continuum. 
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Fig. 5. a: ratio of high flux and low flux spectra (see Sec- 
tion H) to the same power law model. The open circles are 
used for the high flux spectrum and closed triangles for 
the low flux spectrum, b: ratio of the high and low flux 
spectra, c: ratio of high and low flux spectra to power law 
models fit individually to the two spectra, with symbols 
as in panel a. 

6. Discussion 

The standard paradigm for the medium energy X-ray con- 
tinuum emission in radio-quiet AGN is that it is due to 
Compton up-scattering of soft UV photons in some form 
of hot corona (Mushotzky, Done & Pounds 1993 and ref- 
erences therein) . The seed UV photons could be produced 
either by an accretion disk (Haardt, Maraschi & Ghisellini 
1997) or from low temperature (T ~ 3 x 10^ K) thermal 



emission in the corona itself (Torricelli-Ciamponi & Cour- 
In Mrk 766 the broad Fe Ka line seen in the 



1995) 



EPIC spectra provides strong evidence for the presence of 
a cool, optically thick accretion disk (Tanaka et al. 1995). 
This, and the lack of narrow coronal emission lines in the 
long wavelength part of the RGS spectra, favours the ac- 
cretion disk as the source of the UV-soft X-ray bump in 
Mrk 766. 

The UV - soft X-ray accretion disk emission is effec- 
tively bracketed by the soft X-ray excess at ^0.7 keV, and 
the UV emission observed in the CM. The behaviour of the 
(0.4-2.0 keV)/(0.1-0.4keV) hardness ratio (Section|), and 
the X-ray spectral variability (Section ^ implies that the 
soft excess varies significantly less than the Comptonised 
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Fig. 6. a: high flux and b: low flux RGSl spectra, c: dif- 
ference between the high and low flux spectra d: ratio of 
the high flux spectrum to the low flux spectrum. 

continuum. Furthermore, the OM lightcurves are consis- 
tent with zero variability in the UV. It therefore appears 
that the Compton upscattered X-ray continuum varies 
much more than its source of seed photons, the UV - soft 
X-ray bump. This implies that either the Comptonising 
corona is highly susceptible to changes in the soft photon 
field or that the variability of the Comptonised compo- 
nent is intrinsic to the corona itself and is unrelated to 
variations in the soft photon field. Something specific to 
the corona such as changing magnetic fields or the spec- 
trum of particles injected into the corona would therefore 
provide a good explanation for the X-ray variability. 

7. Conclusions 

We have presented RGS, EPIC and OM spectra and 
lightcurves of Mrk 766 from the 60 ks XMM-Newton PV 
phase observation. In all three EPIC cameras the spec- 
tra above 2 keV are well described by a power law with 
a broad Fe Ka emission line. Relative to the power law, 
the EPIC spectra show a deficit of photons between ~ 0.7 
and 2 keV and excess emission at lower energies. The RGS 
spectra contain no strong, narrow emission lines. Warm 
absorber edges are not seen at the expected energies of 
OVII and OVIII in the RGS spectra; an edge-like feature 
at 0.7 keV could only be due to OVII if the absorbing 



material were redshifted more than 10,000 km/s relative 
to the system. Significant variability is detected in the 
X-ray lightcurves and hardness ratios, while no variabil- 
ity is detected in the UV with the OM. A variable warm 
absorber is not required to account for the spectral vari- 
ability, which can be explained by a medium energy con- 
tinuum component which softens as it increases in flux, 
with an additional contribution from a less variable emis- 
sion component at low energy. Changes in magnetic flelds 
threading the accretion disk corona, or the spectrum of 
particles injected into the corona, can explain the observed 
spectral variability. 
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